oo 
o> 

On 



< 



cd 



% 



A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 
2 (12.04.1; 12.07.1; 12.12.1) 



ASTRONOMY 

AND 

ASTROPHYSICS 

1.2.2008 



The effects of source clustering on weak lensing statistics 

F. Bernardeau 

Service de Physique Theorique, C.E. de Saclay, F-91191 Gif-sur-Yvette cedex, France 

February 1, 2008 



Abstract. I investigate the effects of source clustering on the weak lensing statistics, more particularly on the statis- 
tical properties of the local convergence, k, at large angular scales. The Perturbation Theory approach shows that the 
variance is not affected by source clustering at leading order but higher order moments such as the third and fourth 
moments can be. 



1. Introduction 



I compute the magnitude of these effects in case of an Einstein-de Sitter Universe for the angular top-hat filtered 
convergence. In these calculations the so-called Broadhurst and multiple lens coupling effects are neglected. The source 
clustering effect is found to be particularly important when the redshift distribution is broad enough so that remote 
^i | background sources can be significantly lensed by closer concentrations of galaxy sources. The source clustering effects 
^> ■ are shown to remain negligible, for both the skewness and the kurtosis, when the dispersion of the redshift of the 

l/~) | sources is less than about 0.15. 

^H ! 

_ 1 I Key words: Cosmology: Dark Matter, Large-Scale Structures, Gravitational Lensing 

r-: 

Ov 
Or 

q . The construction of gravitational distortion maps, for tracing the large scale structure of the Universe, is a promising 
5h ' tool for cosmology. For the first time, it would indeed provides us with an unbiased representation of the mass 
c/5 distribution in the universe (Blandford et al. 1991, Miralda-Escude 1991, Kaiser 1992). And recent works tend to 
prove that indeed it should be possible to have reliable distortion maps at the level expected for the distortion induced 
by the large-scale structures (Schneider et al. 1997a). 

The cosmological interpretation of such maps is however challenging. The main difficulty is that the background 
sources used to make such measurements are very faint galaxies, whose distances and distribution is not necessarily 
C^ ' well known. It has been pointed out in recent papers (Villumsen 1997, Bernardeau et al. 1997, Jain & Seljak 1997) that 
it was crucial to know with a good accuracy the redshift distribution of those sources. Their distances determine indeed 
the magnitude of the distortion effect: the more distant they are the larger the effect is. The r.m.s. of the distortion 
is however not the only information of cosmological interest to be extracted from distortion maps. In particular the 
departure of the gravitational convergence from a Gaussian statistics in case of Gaussian initial conditions is an 
indicator of the amount of nonlinearity reached by the cosmic density field (Bernardeau et al. 1997, Jain & Seljak 
1997, Schneider et al. 1997b). To get a reliable description of such detailed analysis, one has however to take into 
account the possible effects produced by the intrinsic statistical properties of the sources. 

In this paper I investigate the effect of source clustering for high-order moment of the convergence. In Sect. 2, 
the observational schemes for the distortion and convergence fields are recalled, and the couplings between the lens 
density fluctuations and the background galaxy fluctuations are explicitly written. The calculations are made for an 
Einstein-de Sitter Universe only. In Sect. 3, the implications of such couplings for the skewness, third moment of 
the convergence are investigated. In Sect. 4, the expression of the kurtosis, fourth moment, is computed, taking into 
account the source clustering properties. The numerical applications, in Sect. 5, are made for different models of source 
distributions. 



Send offprint requests to: F. Bernardeau; fbernardeau@cea.fr 



2 Source clustering effects 

2. Observational schemes for the distortion field 

2.1. The filtering schemes 

In the weak lensing regime, and for small-size background objects, the local gravitational distortion effects can be 
entirely described by the local deformation matrix A. This matrix expresses the local linear transform between the 
source and the image plane induced by all the lenses present along a line-of-sight. Its inverse can be related to the 
second order derivative of the gravitational potential through the equation, 



A-\,)=ld-\°- dD^-^^^U (1) 




where Id is the identity matrix, V is the comoving angular distance along the line-of-sight, a(V) is the expansion 
factor and (/^(D, 7) are second order derivatives of the local gravitational potential at the position (T>, 7) along the 
directions i and j orthogonal to the line-of-sight. Expression (|l|) is valid for an Einstein-de Sitter only, but can be 
easily extended to any background geometry. It is written for a given distance of the source, T> Sl that may vary for 
different line-of-sights. 

In practice, the deformation matrix is not directly measurable. The quantities that are directly accessible are the 
shape parameters of the observed galaxies, S 1 , in the image plane. For objects with a large enough extension (compared 
to the width of the point spread function), this matrix can be related to the shape matrix in the source plane, S , 
through, 

S 1 = A- 1 ■ S s ■ A- 1 - }. .. . (2) 

det^- 1 ) v ; 

Then the determination of the direction and amplitude of the ellipticities of the galaxies gives a local estimation of the 
deformation matrixtl. Because of the intrinsic ellipticities of the sources, the cosmic distortion can be detected only 
when a large number of galaxies are taken into account. Averaging over few hundred of galaxies, distortion signals 
down to a few percents are in principle detectable (Blandford et al. 1991). 

The particular quantity that can thus be reconstructed is the local filtered convergence, n. It is directly related to 
the trace of A^ 1 with, 

k = I - tr (A' 1 ) /2. (3) 

This quantity can be obtained directly from the observed galaxy shapes, when it is filtered with a compensated filter 
(i.e. convolved with a function of zero integral, see Kaiser 1995, Schneider et al. 1997b). In general, however, it is 
always possible to obtain a convergence map from a distortion map by solving a differential equation (Kaiser 1995). 
In the following I will therefore focus my analysis on the statistical properties of the local convergence, filtered by a 
top-hat window funtion. 

In the expression (Fy) the relation between the deformation matrix and the gravitational potential is given for a 
unique distance of the source plane. But actually the measured convergences result from averages made over many 
background galaxies that can be at different distances. More specifically the measured local convergence at scale 9q is 
obtained from background objects taken for instance in a solid angle of radius 9q, so that it reads, 



N 
1 

Kg, 



where 7j is the direction of the i th ' source galaxy. The number density of source galaxies that can be used is about 
40 per arcmin 2 (with the usual deep exposures in the I band). For a filtering radius of 20', we expect then to have 
about 50 000 galaxies. This number is large enough to assume that the discretization of the background field is not 
important. We can then work in the coutinuous limit in the source plane. The number density of sources at distance 
D s and in the direction 7 can be written, 

n s (V Sll ) = n s {V s ){l + S s {V sn )), (5) 



1 see however an alternative method proposed by van Waerbeke et al. 1997, called the Pixel Autocorrelation Function, that 
does not use the shape matrices. 
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where n s (T> s ) is the average number density of sourcesa(that fulfill the selection criteria at a given distance), and 
S s (V s ,"f) is their local density contrast. The density n s is normalized to unity, L dT> s n s (T> s ) — 1. Writting the 
equation (|f|) in the continuous limit we have, 

3 /d 2 7 ^ ( 7 ) J*dV s f^ d-D^§^d maSB .(V, 7 )n s (V s )(l + 6 s (V s , 7 )) 

K0 n = — — n ■ (O) 

2 /d2 7 M/ eo ( 7 )/ 2 dP s n s (I? s )[l + ^(2? s , 7 )] 

In this expression, the two filters are a priori both top-hat filters, 

WeM =1 for | 7 |<# ; 

We ( 7 ) =0 for | 7 |>0 O . (7) 

Note however that in general the filters are not necessary the same. For instance one can use a compensated filter 
for the convergence, and still use a top-hat window for the selection of the sources. More generally, instead of giving 
an equal weight to all sources in a given area, it is always possible to weight them proportionally to the inverse of the 
local density. We would then have, 



k 



surf. _ E£L'o w » K (li) 

'o sr^N B ■ 



(8) 



where Wi is a weight associated with each background object. We encounter here a situation similar to the cosmic 
velocity field statistics (see discussion in Bernardeau & van de Weygaert 1996). To get a weight inversely proportional to 
the local density, there are a priori different possible schemes. One can for instance use a two step filtering, with a usual 
filtering scheme on a small grid, and then a final filtering on the larger scale, corresponding to the actual smoothing 
scale of the statistical analysis. One could also think of using similar approaches as the Voronoi and Delaunay methods 
developed by Bernardeau & van de Weygaert. In such a case the weigth wt would be for instance the inverse of the 
surface of the Voronoi cell occupied by a given galaxy. It would lead to a proper 'surface-average' filtering scheme, 



2 J £AV,n,{V,)\l+S,(V„-<)\ 

Unlike in the velocity case, the dependence on the local tracer density contrast does not vanish, because of the 
finite width of the redshift distribution of the sources. In all these formulae both the density contrast of the cosmic 
density field and of the sources are present. They are both random fields and moreover their cross-correlations are a 
priori comparable to their autocorrelation properties. The aim of this study is to investigate the consequence of such 
couplings. 

2.2. The physical effects induced by source clustering 

First of all let me stress that the effects discussed here are different in nature from other couplings that have been 
described previously in the literature. A possible source a coupling between the lens and the source population is, for 
instance the amplification effect that is expected to affect the local number density of detected tracers (Broadhurst et 
al. 1995). This effect is due to the fact that the apparent size of background objects depends on the local gravitational 
amplification making them more, or less, easily detectable. This effect is always present, even if the sources are not 
clustered. It is also expected to vanish when the sources are selected on surface brightness criteria only. Here I assume 
that the sources are, in this sense, perfectly selected, so that the amplification effect is entirely negligible. 

What are then the effects of source clustering? It comes from the fact that the source 'plane' is actually thick and 
inhomogeneous. One mechanism originates from the fact that there might be a significant overlapping between the 
distribution of sources and the distribution of lenses. Let me imagine now for instance there is a large potential well 
at a very large distance, in the overlapping area. Then because the source galaxies trace somehow the matter field, 
you expect to have at the same time a relative excess of close sources. The presence of these sources tends to reduce 
the gravitational signal of the remote lenses. One expects then that the gravitational distortions of the furthest lenses 
will be systematically underestimated. In Fig. 1 a sketch of the actual situation is proposed. 

The implications of such a coupling depends on the level of description one wishes for the local convergence. At 
linear order the expression of the local convergence is not affected. Therefore the variance is not expected to be much 



2 Note that for an Einstein-de Sitter Universe the angular distance at z — > oo equals 2. 
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Fig. 1. Sketch of the lens (thin lines) and source (thick lines) 
density fluctuations. The source distribution is not expected to 
be smooth. Moreover the source density fluctuations are corre- 
lated to the lens density fluctuations. This is particularly impor- 
tant when the overlapping area is large. 



changed at large scale. At small scales, the extra couplings with the source density fluctuations can compete with the 
intrinsic non-linear evolution of the projected density. It should then be addressed with a complete numerical study. 
Another possible effect is due to the fact that the source plane is expected to be 'bumpy', i.e. the average distance 
of the sources may vary from one direction to another. For instance if one observes the gravitational distortion induced 
by a perfectly round potential with a 'bumpy' source plane, the efficiency of the gravitational effect is expected to 
vary from one direction to another creating apparent substructures in the potential well. The importance of this effect 
once again depends on the level of description ones whishes. It is expected in general to create more power at small 
scale, but it can be significant (if it is present at all) only at small angular scale. Perturbatively, this effect is expected 
to play a role only for the kurtosis and moments of higher order. Note that contrary to the previous mechanism, this 
effect is not due to the cross-correlation between the lenses and the sources, but to the intrinsic correlation properties 
of the sources. 



2. 3. Models for lens and source correlations 

The aim of the coming sections is to investigate the implications of source clustering on high-order moments by 
means of Perturbation Theory. Calculations can be pursued only with a model for the mass-galaxy and galaxy-galaxy 
correlation functions. In the following I will assume that a local bias holds and that the local galaxy density contrast 
can be expanded in terms of the linear mass density contrast, 



S a {V,j) =b 1 (V)6£l s (V, 1 )+ b 2 (V) [CLP, 7) 



(10) 



Although there is no complete justification for such an assumption, it is a quite natural from the linear analysis of 
Bardeen et al. (1986) or even from the non-linear descriptionEl proposed by Bernardeau & Schaeffer (1992). Note that 
a skewness of about 3 (Bouchet et al. 1993, Kim & Strauss 1997) for the galaxies implies, 62 ~ 0.5 b\. Furthermore I 
have to assume a given function for the evolution of the bias factors b\ and 62 with redshift. In the following I will 
assume that 



61(D) oc l/a(V) and b 2 (D) oc l/a 2 (V). 



(11) 



According to the results of Bardeen et al. (1986) it corresponds to objects defined with a fixed threshold in units of 
the variance. 

Finally the mass (and galaxy) density field is fruitfully described by its Fourier transform, 



5(V,j) 



d 3 k 

(2tt) 3 / ; 



6(V, k) exp[i(Pk ± -7 + 2? k r )} 



(12) 



where k r is the radial part of the wave vector k and kj_ is its perpendicular part. In the linear regime, the Fourier 
components S(V, k) grow like a(T>) for an Einstein-de Sitter Universe, and they are assumed to obey a Gaussian 
statistics characterized by the power spectrum P(k), 



<5«(©,k)=a(D)5 init .(k); <<W(k)<5i n i t .(k')) 



: (k + k')P(fc). 



(13) 



It implies in particular that the 3-point correlation function of the galaxy field can be exactly factorized in products of 
2-point functions 
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As the numerical calculations will be done at a fixed smoothing scale, it is reasonable to assume that P(k) follows a 
power law behavior, 

P(k) oc k n , (14) 

and at the scales of interest we expect (see Bernardeau ct al. 1997), 

nw-1.5. (15) 

As mentionned before, from a Perturbation Theory point of view, the variance is not affected at leading order 
by the source clustering. In the following, the calculation will concentrate on the skewness and kurtosis of the local 
convergence. 

3. Implications for the skewness 

At large scale, for Gaussian initial conditions, the third moment is given by a combination of the first order and second 
order terms of the local convergence with respect to the initial density field (e.g. Peebles 1980, Fry 1984, Goroff et al. 
1986), 

( K 3 )«3([««] 2 K ( 2 )). (16) 

In the following, it is assumed that the effect of multiple lenses and spurious observational couplings are negligible 
(see Bernardeau et al. 1997). 

3.1. Expressions of the first and second order convergence 

The presence of source clustering does not change the expression of the first order term, 

(D 3 [a* W t \ Pat* f V '^ T> (Ps-T>) <&L(P,7) m , ,,~ 

K «o = 2 J d 7 ^° (7) J dVs J dV V s aJpr ns{Vs) - (17) 

Written in terms of S(k) it reads, 

o r-2 rV 3 p I'D _ j)) r A3]^ 

/ dV s / dV^-^- '-n s (V s ) / -——5 init .(k)exp[ik r V}W(k ± V9o), (18) 

Jo Jo v s J {2-kY' 2 

where the filter W is expressed here in Fourier space. One can further simplify this expression by introducing the 
efficiency function, u)("D), defined by 

u{V) = ^£dV s (Vs v m n s (V s ), (19) 

so that 

4? = -/ dPw(D)/-^i inii (k) cx V [ik r V]W(k ± V9 Q ). (20) 

When the source clustering is neglected, the second order term of the local convergence is given by 

l! -- I &Du{D)l F 2 (k,k') a(V) ft ni t.(k)ft nit ,(k') exp[i (k r + k' r ) V] W[\k ± + k'J Wo], (21) 

where F 2 is an homogeneous function of the wave vectors (e.g. Goroff ct al. 1986). 

When source clustering is taken into account, two extra terms for the second order should be added, 

s.c.(2) (2) 

3fei I' 2 f Vs (V —VXD f d 3 kd 3 k' 

-± j dV s J dVn s (V s ) Ks v ' I S{k)S{k')cxp[i(k r V + k' r V s )}W[\Vk ± +V s k' ± \e } + 

f 2 f 2 /"d 3 kd 3 k' 

h dVu(V) dV s n s (V s ) S(k)S(k')cxp[i(k r V + k' r V s )}W[\Vk ± +V s k' ± \e }, (22) 

In the following, the moments will be calculated with an angular top-hat filter, that in /c-space reads, 

W{k) = 2 J -M, (23) 

where J\ is the spherical Bessel function. 



K B„ - J 
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3. 2. Expression of the variance 

The variance can be calculated straightforwardly from the expression of the linear convergence. Using the small angle 
approximation we have (Bernardeau et al. 1997), 



<- 2 > = 



r[(l-n)/2]r[l + n/2] 



r[l-n/2]r[2-n/2] 



/2] („ +2 ) [ dV „.2 (v]v -(n + 2)- r[(l-n)/2]r[l + n/2] 

-^9 ^dVu{V)V = _______ O 



(n+2) 



(24) 



3. 3. Expression of the skewness 

For a top-hat filter, using the small angle approximation and when the source clustering is neglected, the skewness is 
given by (Bernardeau et al. 1997), 



•S3 



•S3 



(* 3 ) 

(« 2 ) 2 



36 3 (n + 2) 

y 2 



dDw 3 (X>)2?- 2 ( n+2 )o(Z>)//; 



(25) 
(26) 



This result has been obtained from specific properties of the angular top- hat filter (see Bernardeau 1995). For this 
filter no further approximation than the small angle approximation is required. To compute the skewness taking into 
account the source clustering, it is of interest to assume that, 



1 
2W 



sin(0)d0 W (|k + k'|) = W(k) W{k') 



(27) 



where 9 is the angle between the wave vectors k and k'. This property is not exact, but the error it induces is extremely 
weak, less than 1% for n « —1.5. This property implies in particular that the two filtering schemes Eqs. (ph and (H) 
give the same results for top-hat filtering. 

Taking advantage of this expression, we have, 



4°' - S3" 



9 



h dV dV uj{V) 



V(V - V) 
V' 



,.(£>') lu{V) V- {n+2) V'- {n+2) + 



!■&! / dVco(V)n s (V)V-^+^/I 2 , 

L 2 Jo 



(28) 



when source clustring is taken into account. It of course depends only on b\. In Table 1, I give the results for the 
source models described in Sect. 5. 

4. The kurtosis 

In the frame of Perturbation Theory, the expression of the kurtosis is given by 



( K *) c ^( K +)-3(^ = 6( 



,(i) 



,( 2 ) 



) c + 4([ K «l k^) 



(29) 



When the clustering effects is neglected the third order k is given by the integral over the line-of-sight of the local 
third order density. Then, using the general formulae of Bernardeau (1995) one can compute the kurtosis, 



•S4 



Si 



' 2540 
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33(n + 2) + 



21 (n + 2) 2 



dVuj 4 (V)V-^ n+ ^a 2 (V)/L 



(30) 
(31) 



The general expression of the third order moment contains extra terms, 



K s.c.(3) = K (3) _ | / dVf 



T>. 



dV 
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dV s 



dV {Vs V V)V S£ljV)n s (V s )6P(V s ) 
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(32) 



As a result, when one wants to take fully into account the source-clustering effects, many new terms should be 
included. 

In order to simplify the expression of the resulting kurtosis, I introduce the function, 



V(D,V>) = -1h(V - V) {V ' V V)V n s {V) + w{V) n s (V), 

where H is the Hcavisidc function, U{V - V) = if V < V and H{V - T>) = 1 if V > V. 
Then, s| c ' is given by 



(33) 
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Note that, contrary to the skewness case, there is a contributing term 



s bumps_12 6 2 f d p iC , (Pl )p-(«+2) /" dP 2 V(P 1 ,P 2 )2? 2 - ( " +2) / dP 3 V(P 2 ,P 3 )p- ( " +2) c(P 3 ), 
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(34) 



(35) 



which is due to the source auto-corrclation function only. It corresponds to the second mechanism described in Sect. 
2.2. This term does not disappear a priori when the overlapping between the source distribution and the lens efficiency 
function is arbitrarily small. The results shown in table 1 prove however that in practice this contribution is always 
negligible compared to the cross-correlation effects. 

5. Discussions 

5.1. Numerical results 

In order to get numerical results one should choose a specific model for the redshift source distribution. I will assume 
that n* takes the form, 



n s (z) ex / exp [-(z/zo)' 3 ] . 

I make the calculations for three hypothesis, 

model 1 : zq — 1.15, a = 8., f3 = 8.; 
model 2 : zq = 0.75, a = 3., (3 = 1.8; 
model 3 : zq — 0.5, a = 2., f3 = 1. 



(36) 

(37) 
(38) 
(39) 
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Fig. 2. Shapes of the distribution functions of 
the sources (thin lines) and lens efficiency func- 
tions (thick lines) for the three models (models 1, 
2, 3 are respectively plotted with solid, dotted and 
dashed lines). All functions are arbitrarily normal- 
ized to unity 



For the first two models the mean redshift of the sources is about unity, with a larger distribution in the second case. 
The first model has somewhat arbitrary parameters, whereas the second is motivated by results of galaxy evolution 
models (Chariot & Fall, in preparation). It would correspond to galaxies with a / magnitude between 22 and 24. It 
has a redshift dispersion higher than model 1. In model 3, I have assumed a very broad distribution of the redshift 
distribution. The resulting values for the skewness and the kurtosis in such models are given in table 1. As expected 
the first two models give roughly the same answers, because the mean source redshifts are very close. In the third 
model there is a significant number of sources at very high redshift. It lowers the values of S3 and S4. 



Table 1. Skewness and kurtosis with source clustering effects 





model 1 




model 2 






model 3 






(*s) 


1.11 




1.11 






1.50 






A(z s ) 


0.15 




0.42 






0.87 






sr- 


-37.6 + 1.05 61 




-39.9 + 6.4 61 






-29.1 + 7.8 6i 






4' c ' 


2890 - 49 61 + 12 6? - 


-24 62 


3330- 410 61 - 


-217 6? - 


-280 62 


1790 - 390 61 - 


-46?- 


- 290 62 


bumps 

s 4 


1.3 6? 




19 b? 






216? 







Skewness and kurtosis show a similar sensitivity with source clustering effects. For the skewness, the remaining 
corrective term can be as large as 25% of the signal when the redshift distribution of the sources is large. The 
correction is about 30% for the kurtosis assuming that b\ « 1 and 62 ~ 0.5. However for a narrow redshift distribution, 
the corrective terms remain small (about 3% for model 1 for the skewness and 2% for the kurtosis). It shows that 
the high order moments are slightly dependent on source clustering but that this dependence can be controlled if the 
redshift dispersion of the sources is low enough. If one does not want to include pre-knowledge on the galaxy-mass 
cross-correlation to analyze the data, it will be important to reduce as much as possible the width of the redshift 
distribution in the adopted selection criteria. Note that it will be anyway possible to get constraints on the mass- 
source and source-source correlations from counts in cells statistics applied to the selected sources (Schneider 1997, 
van Waerbeke 1997). 

It is also worth to have in mind that these results depend on the shape of the filter. In particular if one does not 
use a top-hat filter the two filtering schemes are not equivalent. This would be the case in particular for compensated 
filters. The "proper" surface-weighting scheme is expected to be, in general, less sensitive to the source clustering. 



5.2. Skewness and kurtosis to measure the non- Gaussian effects 

In the preliminary investigations of the non-Gaussian properties expected to be observable in convergence maps, the 
focus has been mainly put on the skewness. It is indeed the first non-trivial cumulant expected to emerge due of 
mode couplings and its calculation is possible in the frame of Perturbation Theory. However, it would be stupid to 
limit the search of non-Gaussian effects to the skewness only. Even at the level of the shape of the local convergence 
PDF, the skewness cannot entirely characterize the departure from a gaussian distribution. This departure manifests 
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itself by the apparition of a whole set of cumulants. In the context we are interested in, the PDF of k is essentially 
a 2 parameter familly (when a given population of sources is considered) : it depends on the amplitude of the density 
fluctuations and on f2 (neglecting at first view the A dependence) . It implies that all cumulants are somehow connected 
together. In particular the kurtosis S4 appears simply to be an other way of constraining fl (it is in perturbation theory 
independent of as) which means that S3 and S4 are naturally related. This relation can be easily identified from the 
results of Table 1. Indeed one can check that, in the absence of source clustering, 

% * 2.05 to 2.1, (40) 

s 3 

for the three models. The dependence of the ratio S4/s§ on n is expected to be weak, as well as its dependence on 
the cosmological parameters. We encounter here the same situation that has been noticed by Bernardeau (1994) for 
the cumulants of the cosmic density or the cosmic velocity divergence filtered with a 3D top-hat window function. 
This property is a priori valid in the quasilinear regime only. In the strongly non-linear regime the validity of this 
relationship remains an open question although results obtained in a phcnomcnological description of the cumulant 
behavior of the 3D cumulants (Colombi et al. 1997) strongly suggests that such a relation may survive in the non-linear 
regime. The consequences are two fold, 

— S4 might be used as an alternative method of measuring f2; 

— the fact that S4/ 's§ should be approximatly 2 in all cases might also reveal a precious property to use for testing 
the correctness of observational results. For instance, when the bias effects are too strong this ratio tends indeed 
to increase up to about 2.7. 

5.3. Systematics and spurious couplings 

To summarize, the possible spurious couplings that may affect the statistical properties of the convergence maps that 
have been identified so far are: 

— Multiple lens couplings. This is due to the fact that the convergence effects of lenses that are aligned do not 
add linearly. Actually, to compute the effect of combined lenses one should multiply the amplification matrices. 
Departure from the Born approximation appears at the same level of approximation. Both are shown to induce at 
most a few percent effect correction on the skewness (see Bernardeau et al. 1997). 

— Magnification effects. This is due to the fact that the population of selected sources may depend on the local 
magnification (see again Bernardeau et al. 1997). This effect is difficult to quantify a priori. It should be tested 
with the selection algorithms which are actually used. 

— Source density fluctuations. This is one of the effect which is investigated in this paper. It comes from the fact that 
the source plane appears "bumpy" . At large scale, it induces corrective terms at the level of the kurtosis only, but 
it might play a more important role at small scale. 

— Source-lens correlations. This effect, the main effect investigated in this paper, is due to the fact that the sources 
might also somehow trace the foreground lenses when the distribution function of the source redshift is large enough. 
The latter two effects are shown to be negligible when the width of the redshift distribution is small enough (say 
less than 0.15). 

All these effects have been shown to intervene for the high order moments only, in the frame of perturbation theory 
that is at large enough scale (above about 10 arcmin scale). In this approach it is possible to estimate the magnitude 
of these effects and to show that for a reasonable redshift distribution they are negligible. 

In practice however it is likely that the scales of interest will be much smaller than scale at which the perturbation 
theory regime is valid. It would therefore be interesting to extend these results to the small scales. In particular the 
effect of source clustering on the variance might not be totally negligible. 
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